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Electronic transport in highly doped but still insulating silicon at low temperatures is dominated 
by hopping between localized states; it serves as a model system of a disordered solid for which the 
electronic interaction can be investigated. We have studied the frequency-dependent conductivity of 
phosphorus-doped silicon in the THz frequency range (30 GHz to 3 THz) at low temperatures T > 
1.8 K. The crossover in the optical conductivity from a linear to a quadratic frequency dependence as 
predicted by Efros and Shklovskii is observed qualitatively; however, the simple model does not lead 
to a quantitative agreement. Covering a large range of donor concentration, our temperature- and 
frequency-dependent investigations reveal that electronic correlation effects between the localized 
states play an important and complex role at low temperatures. In particular we find a super-linear 
frequency dependence of the conductivity that highlights the influence of the density of states, i.e. 
the Coulomb gap, on the optical conductivity. When approaching the metal-to-insulator transition 
by increasing doping concentration, the dielectric constant and the localization length exhibit critical 
behavior. 



I. INTRODUCTION 

The transition from a semiconductor to a metal upon 
doping is so basic and important that it seems hard to 
believe that this crucial issue is far from being under- 
stood, in spite of the enormous effort invested over the 
years^ 2 . For low concentration, the dopants provide lo- 
calized states which are dispersed in the semiconductor 
crystal. When the host contribution to charge transport 
can be neglected (at low temperature), the system there- 
fore resembles a perfectly disordered solid. At low exci- 
tation energies the charge-carrier transport takes place 
via hopping processes between these spatially randomly 
distributed states and is described by variable range hop- 
ping (VRH). Although the model originally suggested by 
Mott&i neglects correlation effects between the localized 
electronic states, in certain regimes it describes the phe- 
nomena of temperature- and frequency-dependent trans- 
port quite well. Pollak 5 and Srinivasan^ and later in 
particular Efros and Shklovskii (ES) 7 » went one step fur- 
ther and included the long-range electron-electron inter- 
action. During the last thirty years, a large number of 
analytical, numerical, and experimental results have been 
collected^ which clarified the importance of Coulomb 
correlations to some extent. Nevertheless, a complete 
and consistent picture cannot be drawn at this point. 

In a more general view, the disordered localized states 
can be called electron glass: in this terminology Mott's 
model describes the Fermi glass 10 " of non-interacting par- 
ticles, whereas the model taking into account Coulomb 
interactions between the states is commonly referred to 
as the Coulomb glass. 11 With decreasing temperature or 
decreasing electron density, a transition is expected from 
a temperature variation of the resistivity corresponding 
to that suggested for a Fermi glass to a behavior as pre- 
dicted by ES for the Coulomb glass^ The dc conductiv- 
ity is influenced by a Coulomb gap A which opens in 



the density of states around the Fermi edge Ep due to 
the long-range Coulomb interaction. The transition from 
one behavior to the other shows that at higher excitation 
energies the effects of Coulomb interaction are negligible 
and Mott's VRH behavior is recovered. 

For the ac conductivity there exist distinct predictions 
as well. While Mott described the Fermi glass behav- 
ior of non-interacting particles, 4 ES supplemented this 
derivation by introducing the interaction between those 
two states directly involved in the hopping process that 
is triggered by resonant absorption of a photon^ Thus 
with increasing frequency there is also a transition from 
the Coulomb glass regime to the Fermi glass. 

Within the model of ES, 9 the Coulomb gap can be 
roughly evaluated from the temperature dependence of 
the dc conductivity in the two regimes. A more direct 
observation of the Coulomb correlation gap in the den- 
sity of states is possible by tunnelling spectroscopy! 13 i 14 
Quite recently, a controversy arose whether it is really 
the Coulomb gap which determines the crossover energy 
scale from the Fermi- to the Coulomb-glass-like behav- 
ior in the frequency-dependent transport. According to 
ES, this transition is driven by the Coulomb energy of 
the sites forming the resonant pair. However, Lee and 
Stutzmann 1 ^ claimed instead that the transition occurs 
when the photon energy equals twice the Coulomb gap 
width. A subsequent experimental study- 1 ^ and computer 
simulations 1 ^ seem to support the former view. It was re- 
cently suggested 1 ^ that many-particle correlations might 
be crucial for a quantitative description of the hopping 
conductivity. Here we try to resolve this issue by pro- 
viding more accurate data in a wider range of carrier 
concentration and a broader range of frequency. 

A prime example to study the physics of disordered 
solids is phosphorus-doped silicon which can be tuned 
in the insulating regime as well as through the metal- 
insulator transition (MIT) at a critical concentration 
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n c = 3.5 • 10 cm by varying the donor concentra- 
tion n in the single-crystalline host i 1 ' 19 ' 20 When going to 
low enough temperatures and frequencies, the crossover 
from a Fermi glass to a Coulomb glass can nicely be inves- 
tigated in Si:P. 21 Other options to tune the system are 
the application of magnetic fiel d 22 ' 23 and pressur e) 24 ' 25 
as has also been done for SkBi 2 ^ For small enough dop- 
ing concentration n < 2.78 ■ 10 18 cm -3 , a different effect 
was inferred by the strong p(T) dependence; the VRH 
subsides and a simple activated behavior is observed* 1 ^ 
Due to on-site Coulomb repulsion, a Hubbard gap E 2 
splits the impurity band. While a clear distinction be- 
tween uncompensated Si:P and compensated Si:(P,B) is 
seen in dc transport, 27 it is not clear how this effects the 
frequency-dependent transport. 

The work presented in the following mainly concen- 
trates on the temperature and frequency variation of the 
conductivity; the experimental data are compared with 
the theory by Shklovskii and Efros in order to elucidate 
the influence of electronic correlations on the transport in 
disordered electron solids. In addition, the measurements 
on samples with various doping concentrations give in- 
sights into the scaling behavior at the MIT. 



II. THEORETICAL CONSIDERATIONS 

Before we describe our experimental results, we briefly 
present the relevant theories and phenomenological mod- 
els. More details and discussion can be found in well 
established monographs and reviews 4i2i2£ 



A. Fermi glass 

Following Mott's theory of the Fermi glass j* the electri- 
cal dc conductivity cr^c as a function of temperature T for 
a non-interacting three-dimensional electron glass (vari- 
able range hopping) is given by the well-known T -1 / 4 
law 



^(^cxexpl-^Mott/T) 1 / 4 } 

with a characteristic temperature 
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TMott = 21.2 
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Here £ is the localization length and Nq is the density 
of electronic states (DOS) at the Fermi level Ep which 
is assumed to be constant close to Ep in the absence 
of electronic correlations. The Mott temperature Tjyiott 
tends to zero as the MIT is approached Q 

The common theoretical approach to the ac trans- 
port in statistically disordered solids is the so-called 
pair approximation. Each photon-induced hopping pro- 
cess between two localized states has a certain prob- 
ability. The sum over all the hops contributing to 



the absorption process leads to the photon-assisted con- 
ductivity of the system. In the derivation of the ac 
conductivity, a simple one-electron model of the disor- 
dered system is assumed, neglecting any correlation ef- 
fects between the states^ Apart from logarithmic correc- 
tions, the frequency-dependent (w=27r /) conductivity of 
a Fermi glass at T = shows a quadratic behavior: 



o-iH = ae 2 7V 2 ?iw 2 £r 4 



(3) 



Here we denote by o\ the real part of the complex con- 
ductivity a = ai + ioz'-, the imaginary part is related to 
the real part of the dielectric constant ei = 1 — Ana2/u}^- 
The length r u = £ • \n{2Io/hu>} is the most probable dis- 
tance for the hopping processes. Here 1$ is the prefactor 
describing the overlap integral; its value is of the order 
of the binding energy of the electronic stated and hence 
it is commonly taken to be the Rydberg energy of the 
dopant (for Si:P we use Iq ~ 45 meV). The factor a is a 
constant of order of unity. 



B. Coulomb glass 

In the case of appreciable (long-range) electronic inter- 
action, the DOS approaching the Fermi energy is reduced 
due to the Coulomb interaction between the localized 
states, leading to N(E) w (3/n)(E - E F ) 2 e\/e 6 . The 
width of this so-called Coulomb gap is given by 



A ex e 3 7V 1/2 / £l /2 



(4) 



This concept was confirmed by tunnelling experi- 
ment a 13 ' 14 which reveal a parabolic dependence of the 
conductance measuring the DOS. The typical size of the 
Coulomb gap in doped semiconductors is 1 meV. Tak- 
ing the Coulomb gap into account, Shklovskii and Efros 9 
calculated the temperature dependence of the dc conduc- 
tivity, the so-called T~ x / 2 law 



^(^cxexpj-^Es/T) 1 / 2 } 
with the characteristic temperature 
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which is a measure of the characteristic Hartree interac- 
tion strength. At large enough T the accessible energy 
range for hopping processes, A;B(T 3 T Mott ) 1 / 4 , exceeds the 
Coulomb gap A, and Mott's VRH behavior Eq. (fl]) is re- 
covered. This crossover from a Fermi glass to a Coulomb 
glass can be seen as a weak kink in the temperature- 
dependent conductivity^ 

The ac conductivity is again calculated in the pair 
approximation. When taking into account the average 



Coulomb attraction U(r 



e /eir u between two sites 



forming a resonant pair of distance r w , hopping processes 
can even occur between states that are separated ener- 
getically by more than the photon excitation energy two. 
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This long-range Coulomb interaction U{r u ) strongly de- 
pends on the real part of the dielectric constant t\. As 
will be seen later, ei diverges as the MIT is approached; it 
basically screens the Coulomb interaction. The resulting 
optical conductivity 



<7i(w) = ae 2 N% wfrj [hu + U(r u )] 



(T) 



is distinct from Eq. ^ for small excitation energies 
hcu <C U(rJ), when the frequency-dependent conductiv- 
ity yields approximately a linear behavior: 



ae 



a x {w) = NShuC [ln{2VM]" 



(8) 



In the opposite case of weak correlations, hcu ^S> U(r u ), 
the frequency dependence of Eq. |7| is quadratic if the 
logarithmic corrections are neglected. Naturally, it cor- 
responds to Mott's formula Eq. ([3]). Hence the transition 
from linear to quadratic behavior can be considered as a 
transition from probing the interacting Coulomb glass at 
low frequency to the Fermi glass regime at high frequency 
where correlation effects can be neglected. 

Although for both the temperature-dependent dc data 
and the low-temperature ac behavior a crossover is pre- 
dicted, the canonical theory presented so far does not 
represent a direct correspondence. In the former case the 
reduced DOS due to Coulomb interaction - the Coulomb 
gap - modifies Mott's T" 1 / 4 law to the ES T~ x / 2 law. In 
the latter case, however, Eq. (0) and Eq. ||5J) have been 
calculated under the assumption of a constant DOS at 
the Fermi level although meant to describe the interact- 
ing system. Since the DOS of a Coulomb glass exhibits 
the Coulomb gap of width A, this assumption is valid 
only for U{rJ) ^> A, i.e. when the Coulomb interaction 
influences the occupation numbers only to a small extent. 
Then only states outside the gap contribute to ai(u). 

In the opposite case, mainly states within the gap par- 
ticipate in the hopping conduction. For A > U(r u ) > 
huj, Shklovskii and Efrosi^ derive a stronger frequency 
variation than given by Eq. ©. Hence, the incorporation 
of the Coulomb gap results in a super-linear frequency 
dependence 



(TiM = Not , < j /n x 



(9) 



in contrast to the sub-linear dependence of Eq. (jHJl for 
constant DOS. Between both limiting cases a smooth 
crossover is expected. Computer simulations 1 ^ yield 
larger exponents and a more abrupt transition compared 
to the predictions of ESi 9 -^ 

In principle, frequency-dependent measurements pro- 
vide the possibility to investigate electron glasses like Si:P 
in the relevant energy scales. But due to experimental 
difficulties to access the very low energies required (meV 
and less, corresponding to the GHz and THz frequency 
range), the expected crossover from Coulomb glass to 
Fermi glass has been observed only recentl y 15 ' 16 At this 



point it is not clear how the opening of a Hubbard gap in 
the density of states influences the frequency dependent 
conductivity. 



C. Metal-insulator transition 

So far we have only discussed the insulating state of 
disordered systems. In our case the doping of the host 
silicon crystal adds an extra dimension as we approach 
the MIT with increasing doping. This quantum phase 
transition from an insulator with localized states and zero 
dc conductivity at T = to a metal with extended elec- 
tronic states associated with finite dc conductivity can be 
driven by varying external parameters like uniaxial stress, 
magnetic fields or, in our case, the doping concentration 
n i£O£0£i The MIT in joped semiconductors like Si:P 
is due to two facts: the localization arising from disor- 
der (Anderson transition) and from long-range electron- 
electron interaction due to the loss of screening as states 
become localized (Mott- Hubbard transition) ,i i 30 i 32 The 
spatial extent of the electronic states scales with dop- 
ing; the localization length increases as the MIT is ap- 
proached 



£ oc \n c - n\ 



(10) 



where n c is the critical concentration i 9 i 33 i 34 i 35 On the 
metallic side also the conductivity scales; 36 ' 37 <j\ oc 
\n c — n\^. Theoretically, /i is inferred from the correlation- 
length critical exponent \x = (d — 2)v by scaling 
arguments? 3 ^ d indicates the dimensionality of the sys- 
tem. Also the dielectric constant diverges a: 



,37,38 



ei oc \n c 



-C 



(11) 



The power laws are supposed to be universal; that is 
the critical exponents do not depend on the details of 
the geometric structure or the interaction. Previous in- 
vestigations of uncompensated semiconductors inferred 
an exponent /x = 0.5, in contrast to fi = 1 found for 
compensated semiconductors and amorphous metals^ 
However, it turned out that the exponent depends on 
the parameter by which the MIT is tuned and on how 
close to the MIT the experiments are performed. ~ Numer- 
ical simulations and experiments on metal-insulator com- 
posites give values for the critical exponent that scatter 
considerably, for instance £' between 0.5 and I. 39 ' 40 ' 41 ' 42 
In three dimensions, the relation £' = 2/i is expected. 43 
McMillan 44 suggested that both critical exponents are 
related by = v(r) — 1) where 1 < 77 < 3. Applied to our 
case, it is not obvious that the crucial assumption of this 
model really holds that only states within the Coulomb 
gap are of relevance. Recent computer calculations sup- 
port this scepticism^ 

The above considerations neglect the on-site Coulomb 
interaction: For uncompensated semiconductors like 
Si:P, a twofold spin-degenerate ground state Is (Ai) (and 
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higher valley-orbit split states not considered here) is in- 
duced with each P donor atom. Hence uncompensated 
Si:P is always at half-filling of the impurity band; only the 
lower Hubbard band is occupied. With increasing n the 
two Hubbard bands start overlapping at no = 0.871,^21 
but states at the edges are still Anderson localized. It is 
only at n c that they become extended. 



III. EXPERIMENTS AND RESULTS 

The samples for this study were cut from two different 
Czochralski-grown, nominally uncompensated silicon sin- 
gle crystals with a phosphorus gradient along the growth 
axis. 45 From the ingots a number of disks were prepared 
using a diamond wire saw. To remove distorted sur- 
face layers^ the crystals were chemically and mechani- 
cally treated by well established procedures. According 
to Thurber et alM- the doping concentration was deter- 
mined from the room-temperature resistivity employing 
a commercial four probe measurement system (FPP 5000 
by Veeco Instruments). For high doping levels the resis- 
tivity ratio p(4.2 K) / p(300 K) (determined from standard 
dc measurements) is consistent with Ref. |49|. The dopant 
concentration of the crystals used in this study ranges be- 
tween 0.89 • I0 18 cm" 3 and 3.04 • 10 18 cm -3 ; the sample 
properties are summarized in Table [IJ 

For optical experiments specimens of 10 x 10 mm 2 were 
prepared with different thickness ranging from 2 mm to 
0.05 mm and less. The crystals were chemically and 
mechanically polished to optical quality. In combina- 
tion with highly parallel opposite faces, this ensures pro- 
nounced Fabry-Perot resonances. Using a quasi-optical 
Mach-Zehnder interferometer in the THz range of fre- 
quency equipped with backward-wave oscillators as co- 
herent and tunable radiation sources,— the optical trans- 
mission and change in phase could be measured between 
30 GHz and 1.2 THz, corresponding to an energy range 
between 0.1 meV and 5 meV. The Fabry-Perot reso- 
nances due to multi-reflection at the surfaces of the sil- 
icon sample were utilized to enhance the sensitivity and 
accuracy^ From the spectra the real parts of the con- 
ductivity <Ji{uj) and the dielectric constant ei(w) can be 
calculated using the Fresnel formula 2 ^ with an uncer- 
tainty of 10% and 5%, respectively. Complementary op- 
tical experiments were performed using a Fourier trans- 
form spectrometer (modified Bruker IFS 113v) up to fre- 
quencies of 3 THz (corresponding to 12 meV). Depend- 
ing on the dopant concentration, two to four different 
thicknesses were measured of each sample in order to 
optimize the sensitivity. For the sample with phospho- 
rus concentration n — 0.89 • 10 18 cm~ 3 , the absorption 
is too weak to determine the conductivity in the entire 
frequency range. On the other hand, highly conducting 
specimens with n > 2.9 • 10 18 cm" 3 do not transmit suf- 
ficiently to achieve a signal-to-noise ratio large enough 
to analyze the data with respect to <J\ over the whole 
accessible frequency range. 
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FIG. 1: Temperature dependence of the real part of the con- 
ductivity cti for various frequencies and the temperature de- 
pendent dc results for Si:P crystals with (a) n = 1.6-10 18 cm" 3 
and (b) n = 1.97 • 10 18 cm" 3 . 

Temperatures down to 1.8 K were reached with a 
pumped liquid 4 He optical cryostat. As an example, in 
Fig. [2a) the temperature dependence of the conductivity 
as obtained from optical measurements is displayed for 
a sample with n — 1.6 • 10 18 cm" 3 . The ac conductivity 
curves below approximately 10 K suggest that phonon- 
induced charge-carrier transport can be neglected, and 
thus the system is in the so-called zero-phonon regime; 
in other words fiuj > fc^T. This is a basic requirement in 
order to compare the measured conductivity spectra cr(w) 
with the theoretical predictions discussed above. Because 
phonon-induced processes dominate at higher tempera- 
tures, the ac conductivity curves and the dc data coin- 
cide. Obviously, the higher the frequency of a particular 
curve, the higher the temperature at which this occurs. 
Similar behavior is observed for the other samples. Fig. 
QJb) shows the data for a crystal with higher phosphorus 
concentration n — 1.97 • 10 18 cm" 3 . 



IV. ANALYSIS AND DISCUSSION 

A. AC Conductivity 

From the frequency-dependent transmission measure- 
ments at various temperatures, the optical conductivity 
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TABLE I: List of Si:P crystals used in this study. The phosphorus concentration is indicated by n; also the ratio with respect 
to the critical concentration n c — 3.5 • 10 18 cm' 3 is given. p(300 K) is the dc resistivity measured at room temperature; the 
low-temperature to high-temperature resistivity ratio is indicated by p(4.2 K)/p(300 K). The crossover frequency cj c i/27r is 
obtained by the intersection of the linear and quadratic fits of the two regimes. If the frequency-dependent conductivity is fitted 
by power laws eri(a>) oc u> a , the exponents scg and sfg are obtained for Coulomb glass and Fermi glass, respectively. These fits 
give the approximate crossover frequency lu C 2/2tt between both regimes. For samples marked > (<) the conductivity o\ could 
not be determined from Fabry- Perot resonances because the absorption is too high (low). For sample No. 5 only the dielectric 
constant was measured at low temperatures but not the conductivity (denoted by -). To is the temperature of the minimum 
in the dielectric constant, ei corresponds to the dielectric constant measured at T — 1.8 K in the THz range of frequency; \ i s 
the dielectric susceptibility. 
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and dielectric constant could be evaluated for samples 
No. 2 through 6 as described above. This covers the 
doping range from 1.6 ■ 10 18 cm" 3 to 2.57 ■ 10 18 cm' 3 . 
Here we confine ourselves to the lowest-temperature data 
(T = 1.8 K) in order to ascertain that the experiments 
are performed in the zero-phonon regime (ksT < Tiuj) 
for all frequencies of interest. In Fig. [5] the optical 
conductivity for the sample with doping concentration 
n = 1.6 • 10 18 cm -3 is presented. The data can be well 
described by a crossover from a linear to a quadratic be- 
havior of ai(u>) as indicated by the solid lines. As pointed 
out in Sec. [TT1 this corresponds to the limiting cases of a 
Fermi glass with no interaction at high frequencies and 
a Coulomb glass with electronic correlations at low fre- 
quencies, respectively. The corresponding crossover fre- 
quencies w i are listed in Tab. Q] 

The dashed line in Fig. [5] is drawn according to the 
quantitative description given in Eq. ([7]); obviously the 
experimentally observed crossover is much more pro- 
nounced than the smooth transition predicted by the the- 
ory of Shklovskii and Efrosi&2ii2 This strong discrepancy 
of theory and experiment has already been observed by 
previous investigationsj 15 ' 16 

Upon closer inspection as presented in Fig. [3J the low- 
frequency data of all samples are better described by a 
super-linear frequency dependence as suggested by ES 
in Eq. ([9]) than by a purely linear behavior. (Note that 
the logarithmic corrections in Eq. §5§ even lead to a sub- 
linear frequency dependence.) The logarithmic correc- 
tions <7i(u/) oc u>/ \n{2I /hid} for Coulomb glasses have to 
be taken into account if the Coulomb interaction energy is 
smaller than the width of the Coulomb gap, U(r<S) < A. 
The resulting fits are shown for two samples in Fig. [3J 
with the exactly linear fits for comparison. 

In the opposite limit of high frequencies where the pre- 



dictions for a Fermi glass apply, a slightly stronger fre- 
quency variation is recovered than expected from Eq. Q . 
A fit by a power-law a(u) oc lu s yields the exponent s. 
As seen in Tab. HI s(n) increases with doping concentra- 
tion from s = 2.15 for n = 1.6 • 10 18 cm' 3 to s = 2.28 
for n = 2.57 • 10 18 cm' 3 . An overview of the spectra of 
samples with different phosphorus concentration is given 
in Fig. 0]together with the corresponding fit of the super- 
linear behavior at low frequencies and the approximately 
quadratic dependence at higher frequencies. 

We may compare our findings to previous investiga- 
tions deep in the Coulomb glass state. In the radio- 
frequency range (10 2 to 10 5 Hz) at T = 13 mK (in a sam- 



1 




0.001 L-^-J ■ ■ ■ ■ 1 

100 1000 3000 

frequency (GHz) 



FIG. 2: Frequency dependence of a Si:P sample with n = 
1.6 ■ 10 18 cm -3 at T = 1.8 K. The straight lines are fits with 
linear and quadratic behavior respectively. The dashed line 
is a fit with the prediction of Efros and Shklovskii- given in 
Eq. (J7J). The error bars are of the size of the symbols. 
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frequency (GHz) 

FIG. 3: Frequency-dependent conductivity in the Coulomb 
glass regime for the silicon samples with phosphorus concen- 
tration n = 2.29 • 10 18 cm" 3 and n = 2.57 • 10 18 cm" 3 . The 
fits correspond to a super-linear behavior as indicated in the 
graph with Iq = 45 meV. At higher frequencies, the transition 
to quadratic behavior is apparent. Linear behavior is plotted 
as gray lines for comparison. 



pie close to the MIT and under uniaxial stress) Paalanen 
et alw^ found <J\(io) oc lo s with s = 0.9 in Si:P. Castner 
and coworkers investigated Si: As with uj/2it < 2 GHz for 
T > 1 K and found a super-linear frequency dependence 
with 1 < s < 1.5 and a very strong doping dependence 
of s which was interpreted as possible qualitative agree- 
ment with Mott's cu 2 lawj£2 For Si:P they observed an 
increase of s(n) from about unity to 2.2 with increas- 
ing n in a density range comparable to our samples. 53 
Based on previous experiments^ Castner— argues that 
this might be due to a decrease in the Coulomb energy U 
as the MIT is approached. In the light of our results, the 
following interpretation seems plausible: the observed ex- 
ponents for Si:As correspond to the super-linear behavior 
in the Coulomb gap regime, and for Si:P the concentra- 
tion range might include the transition from Coulomb- 
to Fermi-glass behavior but was not identified because of 
the limited number of frequencies in a comparably small 
frequency range. 

As already pointed out above, the transition between 
the two regimes is clearly seen and much more abrupt 
than expected from the general description in Eq. ([7]) 
suggested by ES. From the intersection of the fit curves 
to power-laws, a crossover energy ftw c2 can be defined. 
The values obtained by this procedure are summarized 
in Tab. Q] It should be noted that the resulting crossover 
frequencies lu C 2 are larger than the values lu c i obtained 
from the intersection of linear and quadratic fits (also 
listed in Tab.|J). 

The crossover frequency is of particular interest be- 
cause it is a quantity that can easily be obtained from the 
experimental data and furthermore gives a direct mea- 
sure for the interaction energy U: following Eq. ([7]), at 
the crossover frequency we have Tilo c — U (r u!c ) . In gen- 




frequency (GHz) 

FIG. 4: Frequency-dependent conductivity of the four exper- 
imentally accessible Si:P samples. The respective phosphorus 
concentration is indicated in the graph in units of 10 18 cm -3 . 
The fits correspond to the super-linear and approximately 
quadratic behavior. 



eral one expects the crossover frequency to decrease with 
increasing doping, approaching the MIT. This can be ex- 
plained by the stronger screening of the Coulomb interac- 
tion due to the increase of the dielectric constant, as dis- 
cussed in Sec. lIVBl below. Helgren et alM- suggested that 
the crossover frequency scales with the doping concentra- 
tion ui c oc |1 — n/nd 13 with j3 « 1.65; however, compared 
to those reports we observe a much weaker dependence; 
a power-law fit would yield an exponent (3 « 0.21. Note, 
however, that our values, but in particular the conduc- 
tivity data and subsequently uj c reported in Ref. [H have 
considerably large error bars; hence we should not over- 
interpret this discrepancy. 

Lee and Stutzmann 15 suggested that the width of the 
Coulomb gap affects the crossover energy. Combining 
temperature- and frequency-dependent data on Si:B they 
observe a relation fiu> C 2 ~ 2 A. Taking our sample with 
n = 2.57 • 10 18 cm~ 3 as an example, the temperature- 
dependent resistivity exhibits a kink around 4.6 K which 
indicates the transition from the Mott regime to the 
Efros-Shklovskii regime. From the fits in both limits 
[Eqs. (P) and ©] we obtain T Mo tt = 2.25 • 10 6 K 
and Tes = 790 K which yields a Coulomb gap A = 
1.15 meV, following Ref. [LJ. The corresponding fre- 
quency 280 GHz is around half of the crossover frequency 
iv c /2n — 570 GHz, as suggested in Ref. [H. 

Whether this agreement indicates a close link between 
the width of the Coulomb gap and the crossover between 
the two regimes indicated in ac transport, however, re- 
mains an open question. One reason is that the validity 
of Eq. ([7]) is unclear for the important and experimen- 
tally relevant case of U(r u ) k 2 A. Lee and Stutzmann 
proposed that the sharp crossover is controlled not by 
the average interaction strength U as in Eq. ((TJ, but in- 
stead by the Coulomb gap. In other words, it is not the 
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mean Hartree energy between sites forming a resonant 
pair that is relevant, but the correlation energy 2 A. For 
this reason the single-particle gap measured by tunnelling 
is much larger than the renormalizcd Coulomb gap that 
governs transport. Helgrcn et al., 16 on the other hand, 
suggest that the Coulomb interaction determines the ob- 
served frequency-dependent crossover from ES- to Mott- 
like hopping conduction occurs. Based on recent com- 
puter simulations, Basylko et alM also argue that the 
transition is driven by the Coulomb energy of sites form- 
ing resonant pairs and not by the width of the Coulomb 
gap. They calculated the number of sites outside the 
Coulomb gap relative to the total number of sites par- 
ticipating in the ac conductivity at a given frequency 
u), and found that within the frequency range of inter- 
est only sites inside the Coulomb gap contribute to the 
ac conductance, in accordance with our observation of a 
super-linear frequency dependence. 




n(10 18 cm" 3 ) 



FIG. 6: Dielectric constant ei of Si:P as a function of phos- 
phorus concentration measured at T = 1.8 K. The solid line 
represents a fit by ei — ehost oc (n c /n — l) - ^ with f" = 1.2. 



B. Metal-Insulator Transition 

Besides the frequency-dependent conductivity, the di- 
electric constant ex(u>,T) was measured for each Si crys- 
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FIG. 5: Temperature dependence of the dielectric constant ei 
for Si:P samples with different phosphorus concentrations n 
as indicated. The curves correspond to the left and right axes 
as indicated; also note the different temperature axes for the 
two frames. The data points are mean values averaged over 
the frequency range investigated. 



tal of different phosphorus concentration. In the entire 
frequency range 30 GHz < u/2ir < 3 THz, the dielectric 
constant £i(w) is basically frequency independent up to 
about T = 30 K. 

Surprisingly, as presented in Fig. [5j we find a non- 
monotonic temperature dependence e±(T) with a pro- 
nounced minimum at low temperatures To which might 
be an inherent property of the zero-phonon regime.— 
Depending on the phosphorus concentration, the min- 
imum decreases from approximately Tq = 11 K for 
n = 0.89 • 10 18 cm" 3 to T = 5 K for samples with 

-3 



n = 2.91 • 10 cm -3 . Within our range the concentration 
dependence of To can be approximated by a linear behav- 
ior. Temperature-dependent measurements on undoped 
silicon do not exhibit such a minimum in the dielectric 
constant.— Also for samples with higher concentration 
there is no minimum evident in the observed tempera- 
ture range; for these samples To is very likely below the 
lowest temperature of T = 1.8 K that could be reached 
with our optical cryostat. It is tempting to suggest that 
the opening of the Hubbard gap E^ leads to an additional 
contribution to the dielectric constant Aei oc E^ 2 ■ Sim- 
ilar experiments on compensated samples could clarify 
this issued 

As can be seen in Fig. for T > To the dielectric 
constant increases with temperature because relaxation 
processes due to phonons become more likely. Note that 
the values of the dielectric constant for each temperature 
are obtained as the average of dozens of data points ac- 
quired in the frequency regime under inspection; they are 
very accurate (better than 0.5 %) and geometrical effects 
can be ruled out. 

The low-temperature (T <C Tq) dielectric constant e\ 
diverges when the phosphorus concentration increases to- 
wards the MIT as depicted in Fig. [51 In the present 
case it is more appropriate to subtract the lattice term 
£host = H.7 of undoped silicon in order to gain direct in- 
sight to the critical behavior of the disordered localized 
system itself. Consequently Fig. [7] shows the concentra- 
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FIG. 7: Dielectric susceptibility An\ of Si:P crystals versus 
the reduced doping concentration 1 — n/n c with n c = 3.5 ■ 
10 18 cm -3 at T = 1.8 K. 



tion dependence of the dielectric susceptibility fitted by 

47rx = ei - e ho st oc (1 - n/n c y c . (12) 

The critical exponent £ = 1.68 is in reasonable accord 
with previous reports (C = 2.0, Ref. LL6J) taking the er- 
ror bars of about 20% into account when comparing dif- 
ferent samples. Significant deviations are observed for 
small phosphorus concentrations, i.e. further away from 
the MIT. 

A resonable fit is obtained in the range 0.13 > 1 — 
n/n c > 0.5. It should be pointed out that there is no 
clear indication of the opening of the Hubbard gap that 
was inferred to occur around 1 — n/n c 0.8 from trans- 
port measurements i2I We also note that for more heavily 
doped samples on the metallic side, the critical region of 
the critical exponent /i m 1 of electrical dc conductivity, 
<Jdc <x| n — n c | M , is much smaller, i.e. 1 — n/n c ss 0.2. 19 

There is some confusion in the literature because of 
an apparent "scaling" with the variable n c /n — 1, i.e. 

X oc \n c /n — ll - '' , extending to much larger reduced 
concentrations . 54 ' 58 However, there is no theoretical jus- 
tification for such a "scaling" . For demonstration pur- 
pose only we note that the line shown in Fig. [5] is actu- 
ally a fit with this function, yielding £ = 1.2 in strong 
contrast to £ = 1.68, obtained from the fit to Eq. (fT2"j) . 

From temperature-dependent dc measurements Hor- 
nung et al£^ evaluated the concentration dependence of 
the Mott temperature Tjyiott and via Eq. ([2]) the scaling 
exponent v = 1.1 for the correlation length. A scaling be- 
havior is expected only very close to the MIT, neverthe- 
less we can extract the relative concentration dependence 
of the localization length from the frequency-independent 
prefactor in Eq. ([5]), if we take the low-temperature value 
of the conductivity and dielectric constant. These results, 
plotted in Fig. [5J can be well described by a critical ex- 
ponent of v = 0.87 according to Eq. (fTTJ)) . Our findings 
are in excellent agreement with similar experiments for 
which v = 0.83 was reported^ They fulfill the Harris 
criterion 5 ? v > 2/d, with d the spatial dimension, which 




reduced concentration 1-n/n 

c 



FIG. 8: Scaling behavior of the localization length in Si:P. 
The critical concentration is n c — 3.5 • 10 18 cm -3 . 

implies that the disorder does not affect the critical be- 
havior. 

Castner and coworker s 28 ' 53 suggested a relation of the 
dielectric susceptibility to DOS and to the localization 
length 

X = e 2 N(E F ) i 2 . (13) 

Thus the ratio of both critical exponents (,/v is expected 
to be 2 which is perfectly confirmed by our measure- 
ments: C,/v = 1.97. Looking at related systems, like 
doped germanium, indicates that the observed behavior 
is a general one, but also that the actual exponents de- 
pend on the doping range under consideration and the 
degree of compensationj 60 ' 61 

V. CONCLUSIONS 

Our investigations of the electrodynamic properties of 
phosphorus doped silicon in a broad frequency range from 
microwaves up to the far infrared underline that elec- 
tronic correlations between the localized states play an 
important role at low temperatures. The crossover in 
the optical conductivity from a linear to a quadratic be- 
havior as predicted by Efros and Shklovskii is observed 
qualitatively; however, the simple model does not lead to 
a quantitative agreement. Our measurements show that 
the crossover energy Tux> c is not solely determined by the 
interaction energy U(r u ). Furthermore, the Coulomb gap 
in the density of electronic states explains the observation 
of a super-linear frequency dependence in the Coulomb 
glass regime. When approaching the MIT with increas- 
ing doping concentration, the dielectric constant and the 
localization length show critical behavior, in good agree- 
ment with theoretical predictions. 

However, there remain open questions to be addressed 
in the future. Concerning the dielectric constant, there 
is the unclear origin of the minimum in the temperature 
dependence. As far as the conductivity is concerned, our 
results should be confirmed in comparable systems, like 
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Si:B or Si:As. Of superior importance is a detailed in- 
vestigation of the influence of compensation: while the 
theory of Shklovskii and Efros was done for compensated 
semiconductors, most experiments (ours as well) have 
been performed on nominally uncompensated samples. 
Furthermore, studies in even larger parameter ranges are 
highly desired to close the gap not only in frequency but 
consequently also in concentration between the GHz re- 
sults of Lee and Stutzmani*i£ and the THz results of the 
present and other recent studies^ Here the concentra- 
tion dependence of the frequency crossover is of particu- 
lar interest. 

Finally, it would be of interest to analyze the 
frequency-dependent conductivity at elevated temper- 
ature, i.e. when the zero-phonon regime is left. For 
ksT m hu> both photon-assisted and phonon-assisted 



processes are equally important. Although the data ex- 
ist, suitable models are still lacking. 
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